Thermo-optic noise in coated mirrors for high-precision optical measurements 
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Thermal fluctuations in the coatings used to make high-reflectors are becoming significant noise 
sources in precision optical measurements and are particularly relevant to advanced gravitational 
wave detectors. There are two recognized sources of coating thermal noise, mechanical loss and ther- 
mal dissipation. Thermal dissipation causes thermal fluctuations in the coating which produce noise 
via the thermo-elastic and thermo-refractive mechanisms. We treat these mechanisms coherently, 
give a correction for finite coating thickness, and evaluate the implications for Advanced LIGO. 
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I. INTRODUCTION AND MAIN RESULT 

^^ I Thermal fluctuations in the coatings used to make high-reflectors are becoming significant noise sources in precision 

~j ■ optical measurements [1][2][3]. Though masked by other noise sources in the currently operating first generation 
] interferometric gravitational-wave antennae (e.g., GEO [htt], LIGO [htt], TAMA [htt], Virgo [htt]), designers of 
^ . second generation gravitational-wave antennae expect coating thermal noise to be the dominant noise source in the 
K~5 ' detector's most sensitive frequency band [4]. Reduction of coating thermal noises has the potential to significantly 
f— ^ , increase the sensitivity, and thus the detection rate, of these large scale detectors. 

^f) '- Coating thermal noises are defined by differences between the coating material and the substrate material[17]. 

There are two recognized sources of coating thermal noise, mechanical loss and thermal dissipation. The first of these 

' - -j'' leads to "coating Brownian" noise, which, while not the topic of this paper, serves as measure against which we will 

^^[ compare our results [5]. The second, thermal dissipation in the coating, leads to temperature fluctuations, which 

I . can cause "thermo-optic" noise via thermal expansion of the coating, and thermal change in refractive index of the 

{^JT)' coating material [6]. 

Despite their common origin, coating thermo-elastic and thermo-refractive noises have not been treated in a coherent 
way [7] [8] . Since the two mechanisms can be of the same order of magnitude a coherent treatment has the potential 
to greatly change the predicted magnitude of thermo-optic noise. 
^L ' The purpose of this paper is to unify the thermo-optic mechanisms. The formulaic result of this unification is 

r~-«, I presented later in this section, and derived in section IL A correction for coatings of non-negligible thickness is given 
t^^ ■ in III. In section IV, we evaluate the thermo-optic and Brownian noises expected to be present in Advanced LIGO, 
^^ ' given current understanding of coating material parameters and detector design. Finally, in the appendices we give 
r — . equations for evaluating the average material constants of a multi-layer coating, we describe the dependence of the 
f^ ' reflection phase of a coating on its temperature, and we relate coating thermo-elastic noise to substrate thermo-elastic 
(>D , noise. 

^D • The power spectrum of thermal fluctuations responsible for thermo-optic noise, as observed by a sensing beam with 

a Gaussian proflle, is given by[8] 



> 
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(see the table at the end of this section for a list of symbols representing material parameters, their deflnitions and 
units) [18]. 

These thermal fluctuations result in fluctuations in the phase of a field reflected by a mirror's coating, which for 
the sensing beam arc equivalent to changes in the inferred position of that mirror via the simple relation[19] 

dT ~ d<l> df~ 4^aT ~ 'f'df' ^' 

where we define di = —\/Att to avoid repetition. 

The spectral density of thermo-optic noise in a mirror's measured position is given by 

qAz qAT I az^^'^c - ,C!c\ ,^ 

Sto = Sto ( "0"^ ~ ' CT ) ^ ' 

which accounts for both thermo-optic mechanisms in d(j)c/dT, the coating's overall reflection phase sensitivity to 
temperature [20]. Thermo-optic noise is explicitly limited to the coating by taking the difference between the temper- 
ature sensitivity of the coating and that of the substrate {d<j>s/dT = Sia^, integrated over the coating thickness d, 



and weighted by the relative heat capacity to give the term in equation 3). Here, we will continue without diversion 
to an expression which can be easily evaluated, but to get a more precise result the thickness of the coating must be 
corrected for as described in section III, and the value of d(j)c/dT computed as described in appendix B. 

Thermo-elastic expansion of the coating is complicated by the mechanical constraint of its attachment to the 
substrate [9] . Under the assumption that the coating elastic coefficients are similar to those of the substrate material, 
the effective thermal expansion coefficient simplifies to 



ttc -- 2ac(l + crc) 



(4) 



The same expression is valid for the constrained thermal expansion of the bulk material in a semi-infinite substrate, 
cis- A more complicated expression for cuc applicable in the case of differing elastic coefficients is given in appendix A. 
Both dn/dT and thermal expansion play a role in changing the reflection phase of a coating, so while this mechanism 
is refereed to as "thermo-refractive" , this is something of a misnomer. Nonetheless, for a high-reflection coating made 
of 1/4-wave doublets with a 1/2- wave cap layer (see figure 1), the thermo-refractive mechanism can be expressed in 
terms of the equivalent dn/dT of a hypothetical single layer of n = 1, a = material backed by a perfect reflector. 
In this case, with a hypothetical layer of thickness A, we find 
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where Bx is the fractional change in optical path length with respect to temperature in material X 

Bx ^ Px + axnx 
with X G {L, H} either the low-index material L, or the high-index material H. 

y2-wave cap %-wave doublets 

\ ^ ^ 
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FIG. 1: A high reflection coating made of 1/4-wavc doublets with a 1/2-wave cap layer. 

Accepting the approximations above, we can combine thermo-elastic and thermo-refractive (TE and TR) mecha- 
nisms to write the total coating reflection phase sensitivity to temperature as 



which allows us to rewrite (3) as 






Sj-Q ~ SrpQ ( acd — /3A — agd-p^ 



Cs 



(7) 



(8) 



Though this paper includes some refinements to previous works, the primary result is that the relative sign between 
the TE and TR mechanisms is negative. 



A. Symbol Definitions 

The physical constants, material parameters, and frequently used symbols in this paper are: 



symbol 


I name 


SI unit 


kB 


Boltzmami's constant 


J/K 


T 


mean temperature 


K 


LU 


angular frequency 


rad/s 


C 


heat capacity per volume 


J/Km^ 


K, 


thermal conductivity 


W/mK 


n 


refractive index 




a 


thermal expansion 


1/K 


P 


dn/dT 




E 


Young's Modulus 


N/m^ 


a 


Poisson ratio 




A 


beam wavelength 


m 


ro 


beam radius (1/e^ power) 


m 


d 


coating thickness 


m 


dl 


d I^z/d(j) = -A/47r 


m 



Material parameters that appear with a subscript refer to either the bulk substrate material parameter, subscript 
s, the average coating parameter, subscript c, or to one of the coating materials L, for low-refractive index, or H, for 
high-refractive index. 

Material parameters which appears without a subscript, but as a function of r take on the value of the material at 
the location r. Thus, a{'r) is cic when r describes a point in the coating, and a^ for points in the in the substrate. 

Bars are used above symbols to express an "effective" coefficient. These coefficients have the same units as their 
bar-less counterparts and the same general meaning, though taken in a specific context. For example, a has the same 
units as a, and is a thermal expansion coefficient, but only in the context of a semi-infinite medium. 

II. REFLECTION PHASE NOISE 

This section will derive equation 3 from the Fluctuation-Dissipation Theorem (FDT)[10]. We'll start by performing 
the derivation of Fejer's result for thermo-elastic noise using Levin's simpler approach [9] [8]. The solution to the more 
general problem of thermo-optic noise is derived second, following the same path. 

The thermal fluctuations which are the source of thermo-optic noise are important to optical measurements because 
they change the result of position measurements based on reflecting a field from a mirror. The fields used in these 
measurements are well described by a normalized Gaussian intensity profile [21] 



/(rj_) =^e-2'-xAo, (9) 

where tq is the beam radius, and r^ = a;^ + y^ is the radius perpendicular to the beam's propagation direction (along 
the z axis). 

To go from thermal fluctuations to measured displacement noise we return to the foundation of this analysis. Our 
application of the FDT starts with a gedanken experiment in which we consider an oscillating power injection in a 
small volume SV located at r 

JL=TFosmiLot)qir}. (10) 

where Fq is an arbitrary scale factor, and lo is the frequency of interest [22]. The form factor q(f) connects the 
measurement variable z to temperature fluctuations 5T{r^ t) in the mirror via 



/ 



5T{r,t)q{r) (11) 

where the integral is formally over all space, though the integrand is presumably zero outside the mirror and its 
coating. 

Power injection leads to heat flow and thus dissipation as expressed by 



W^{^ -{V5TY), (12) 



where the average (...) is over cycles of the power injection. FinaUy, the FDT relates this dissipation to the spectral 
density of noise in the associated measurement variable by 



,A, _ 8kBTW 



S^^ = -^^;^. (13) 





In the next section, as a illustrative example, we will derive Fejer's result for thermo-elastic noise using the approach 
outlined above. The same approach is applied to the more complicated problem of thermo-optic noise in section II B. 

A. An Example: Coating Thermo-elastic Noise 

As a concrete example, we will first apply the above formalism to derive coating thermo-elastic noise in the absence 
of any thermo-refractivc mechanism, (previously performed in [7] [9] ) . The thermo-elastic readout variable 

ZTE = <{ ST{r,t)qTEir) (14) 

describes the sensing beam's averaging of the thermally induced displacement of points on the mirror's surface. The 
solution to the non-trivial problem of thermal expansion of a thin coating on a semi-infinite substrate is presented in 
[9], appears in our equation A2, and is contained in the effective thermal expansion coefficient a. From this we can 
write simply 

qTsir) = I{r±)a{r), (15) 

which leads to a thermo-elastic power injection 

Pte 



6V 



TFQsm{ujt)I{rj_)a{f}. (16) 



To remove the component of power injection which results in little temperature gradient and thus little heat flow, 
we subtract the substrate contribution 

1 Pte^ 1 Pte 1 Pte^ 



C{r) SV C{f) 6V Cs 6V 

= Tf...„M,7(.j(||-g 

We can then recast this into the form of (10) as 

Pte. 



^^ TFosm{u;t)qTEAr), (17) 

where we have identified the coating thermo-elastic readout form factor 

QTEArO - Hr±) (a{r) - a.^) (18) 

which is zero in the substrate by design[23]. 

To maintain the simplicity of this example, we will assume that the coating and substrate are uniform, and that 
the coating is of thickness d which is small with respect to the thermal diffusion length 

With this assumption, we can consider all energy to be generated in this thin layer at the surface of the substrate and 
then flow inward. Integrating (17) over z we compute the energy flux into the substrate to be 

^^^ = TFa sin(wi) /(r_L ) Aad (20) 

oA 



where 



C 
Aa = ac-as-^. (21) 



In order to connect this heat injection to W in (12) we solve the diffusion equation 



C^ = kV^ST (22) 



with the boundary condition that the injected energy flows inward 

(23) 



Pte dST 



SA " dz 



z=0 



which ignores the very small radiation loss, as in [8]. 

Further assuming that re ^ rx, we can ignore diffusion in the transverse dimensions, which yields the solution 



dST -TFn 



z 



^ _ ig^-T sin u:t -= — /(r^) Aad, (24) 

dz Ks \ V2rT, 



from which we can compute the power dissipation 

Ks (d6T\ 



W- 



TEa 



T \ dz J 



Finally, returning to (13), we arrive at the coating thermo-clastic noise spectrum 

A. 2V2kBT^ ^^_^^2 

^TE, = 2 / ^ (^"O) 

= S^Si^adf (26) 

which is equal to that of [9], and [7] under their simplifying assumptions. 

B. Coating Thermo-optic Noise 

The thermo-elastic noise described above assumes that the relevant readout variable is based on the position of the 
surface of the mirror. Interferometric sensors are, however, actually sensitive to the reflection phase of a surface as 
well as its position. 

In the case of reflection from a planar surface, the position and reflection phase are related simply by 5z ~ 
—6(j) A/47r — 9? (50, but for multi-layer coatings the relationship can be more complicated. To account for this we 
generalize (14) to yield the thermo-optic readout variable 






(27) 



where as before we have subtracted the substrate contribution so as to remove the component of heat injection which 
can be handled adiabatically. From this we identify the thermo-optic form factor 



which is, as before, zero in the substrate. 
Plugging into (10), we get 



''-^-TF,sH^t)Iir.)dl('^^'^ ^(^)' 



5V ^ ' -^ ^' f y Qrp Cs dT 



Again we assume that the coating and substrate are uniform, and d ^ r^ <C rg, so we can integrate over z to get the 
energy flux 

?^ = TFo sin(^t) /(r^) [df-^ - a.dg) (29) 

where 

^ = /i Mf) (30) 

dT Jo dT ^ ' 

is the overall reflection phase sensitivity of the coating to temperature, as described in appendix B. 
Following the path used for (26) above, we arrive at 



and thus 



which matches (3). 



S^6-S^S{d;^-asd^X (32) 



III. THICK COATING CORRECTION 

Here we will allow for finite thickness coatings by removing the assumption that d ^ r^, while continuing to assume 
re 2> rx [24] . To do this we will need to solve the heat diffusion equation accounting for power deposition and diffusion 
in the coating. Generalizing (22) to include a source term, but limiting heat fiow to the z axis 

dST dHT P 



with the one-dimensional power injection 

Pro _ 1 Pro 
~5L ~ I{r±) ~SV~ 



(34) 
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dT "C, 



We will approximate the thermo-optic power deposition in the coating with a constant thermo-elastic component, 
and a Dirac delta function for the thermo-refractive component since its effect is limited to the first few layers of the 
coating. We can express this as 

^^° - TFo sm{ujt) (a, - 5{z) Ptr - a^^) (35) 



oL \ L: 

where we define 

PTR^a^d^dl^^pX (36) 

with d4>c/dT and /? as given in appendix B. 

Following the method used in [9], we transform (33) to a second-order differential equation in z 



where the relationships between the new and old variables are 



dT{z,t) ^ ne{e"^'e{z)) (38) 

P{z,t) 



6L 



= LuCne{-ie''^'p{z)) 



.LOC \/l 

The homogeneous solutions to this equation in the coating and substrate are 

6hc{z) = 6'dc cosh (7c z) 
Bhs{z) = 6*^8 exp (-7^2;) 

where the coefficients dd will be determined by boundary conditions at z = d. These equations satisfy the boundary 
conditions of no heat flow at z = or z = 00 . 

The particular solutions needed are for the two kinds of sources, TE and TR, both of which are limited to the 
coating. The source terms are 

PTE = —TT [(^c-as— ] = — ^Aa 

1^ F 
Ptr{z) = -^{z) —pT- Ptr = 5{z) PTR 

and the corresponding particular solutions are 

SpTE = —PTE 
0pTR{z) = ~Jc PTRCXp{-'-fcZ). 

We put all this together with boundary conditions at z — d that ensure continuity of temperature and conservation 
of energy 

Ohs{d) = 0pTE+OpTRid)+0hc{d) 

d d 

KsTrdhs{d) = Kc^r-{9pTR{d)+ehc{d)) 

oz oz 

Ode = Pte + PTfl7cexp(-7cd)(l - R)/'ipd 

0ds = -i?exp(7<;d) {pTE sinh(7cd) + ptrIc) l^d 

ipd = cosh(7crf) + i?sinh(7cd) 

„ KrCr Hrlr HrfTt 

R = _L^^_m ^_lJ^_ (39) 

V KsCs Ksjs KsTTc 

Before we lose ourselves among the equations, recall that our goal is to find the time averaged dissipation W, which 
is related to the temperature gradient in equation 12. We now have 6{z) in hand, and (38) relates this to ST, so our 
destination is near. Summing the homogeneous and particular solutions to get 0{z), and taking the derivative with 
respect to z, wc find 

d6c(z) 

- 7c(6'dcSinh(7cz) + 7c PTi?exp(-7cz)) 



to find 



dz 

dz 

From equations 12 and 38 we can see that 



-^sOds exp (-7s2) 



Txrthick 

WrpQ 



1 / 


n ^ 


fdST 


V 


nrl\ 


V't{3z) 


1 


I?'r 


deiz) 

dz 


2 


2nrl 





(40) 



where the transverse integrals over I{r±) ^ have aheady been performed. 

To arrive at a correction factor for thick coatings, we normahze the corrected thcrmo-optic dissipation above by 
that of a thin coating given in equation 31, 

TJ/thick Q^Zthick 



Wto S^, 



TO 

Taking the integral over the coating and substrate, we end with a complicated expression for the correction factor 

To = 2(sinh(0-sin(C)) + 2i?(cosh(0-cos(e)) 

Fi = 8sin(^/2)(Ecosh(^/2)+sinh(^/2)) 

Fa = (l + i?2)sinh(C) + (l-i?^)sin(e)+2i?cosh(0 

Td = (l + i?^)cosh(0 + (l-i?2)cos(0 + 2Esinh(0 



where we have made the following substitutions 



Ptr. dpTE ,.„^ 

VR = -, -^ — , PE = -, — — (43) 

apTE + PTR. apTE + PTR 



using the dimensionless, frequency dependent, scale-factor 



C = ^ = ,l^^d. (44) 

TTc V Kc 

Note that the power deposition fractions pe and p^ can also be written as 

~0X Aad ,^^, 

Applying this correction to equation 8 gives 

Sto = S^S Tic (Aad - ^A) ' . (46) 

For d <^ rxc or ^ <C 1 , we can use the much simpler expansion 

3i? ^ 

which goes to 1 as ^ goes to 0. In the case of a very thick coating, with d ^ rxc, the thermal fluctuations which 
generate noise via TE and TR mechanisms become independent, and thus they add in quadrature [2 5] 

T^ ^ 2p| , pj, (^g) 



Rii+R)e R 

Thus, this correction expands our understanding beyond the simple notion that the TE and TR mechanisms have a 
relative negative sign. Now we can say that TE and TR mechanisms have a relative negative sign if d ^ ttci are 
partially coherent and partially canceling if d ~ I'Ta and act as independent noises if d 3> ttc- 

IV. IMPLICATIONS FOR ADVANCED LIGO 

Having clarified the relationship between the thermo-optic mechanisms, a recomputation of the impact of this noise 
source is in order. We will also take this opportunity to use the most recent information about the physical properties 
of the materials involved, and to apply an additional correction factor for the less-than-infinite size of the mirror. To 
highlight the implications of this work, the results will be compared with Harry's result for coating Brownian noise 

[5]. 
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FIG. 2: Thermo-optic noises and Brownian noise for an Advanced LIGO input mirror. 
ETM Coating Noise: N doublet = 19, T = 5.0 ppm 
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FIG. 3: Thermo-optic noises and Brownian noise for an Advanced LIGO end mirror. 



The Advanced LIGO mirrors are high reflectors with a multi-layer coating of alternating Si02 and Ta205. The 
input mirrors will have a power transmission of T = 1.4% with tq = 5.5 cm, while the end mirrors will have T ~ 5 ppm 
with ro = 6.2 cm. The mirrors are made of fused-silica, are 34 cm in diameter and 20 cm thick for a total mass of 
40 kg. In figures 2 and 3 we plot the coating related noises for coatings made of 1/4- wave doublets. 

We take the finite test- mass correction from [11] which, with the mirror and beam-size parameters given above, is 
Cfsm — 0.98. This multiplicative factor affects only the thermo-elastic mechanism, 



Aa 



fsm 



Cfsrn^a, 



(49) 



as it represents a bending of the optic due to strains produced by the coating. Adding this correction to equation 46 
gives 



-'TO 



Sxo Ttc {Aafsmd - (3X) . 



(50) 



where we use Aafsm in (45) when computing Ttc- 

Figure 4 shows a representative Advanced LIGO sensitivity curve [26]. While the difference between the result of 
equation 50 and a conservative estimate which simply takes the sum of the TR and TE mechanisms is less than 10%, 
our coherent treatment of TO noise makes clear that it should not be considered a driving force in Advanced LIGO 
coating design. 
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Advanced LIGO Sensitivity 




FIG. 4: An Advanced LIGO sensitivity curve. The themio-optic curve labeled "NEW" uses equation 50, while the "OLD" 
curve uses a conservative estimate of TO noise: the sum of TR and TE, with the TE correction factor of 1.56 from [11]. 

It should also be noted that some of the material parameters used to make these figures are poorly constrained. 
The thermal conductivity of Ta205 is simply assumed to match that of sapphire [9]. Fortunately this only effects the 
thick coating correction factor, and no reasonable value significantly changes the result below 1 kHz. The value of 
/? for Ta205 is also poorly constrained, but again the range of tolerable values is large. Thermo-optic noise remains 
below the conservative "OLD" curve in figure 4 for values between — 10^"'/K and 3 x 10^"'/ K. The value of /Jr^^^ q^ 
used herein, from [Gretarsson] , is comparable to previous values [12] [13]. 

The values of material parameters used for figures 2, 3 and 4 are: 

symbol Ta205 unit 



a 

K 

c 

E 
a 



3.6 

14 

33 

2.1 

140 

0.23 

2.06 



10-«/K 

10-6/ K 

W/mK 

GJ/Km^ 

GPa 



symbol Si02 unit 



a 

P 

K 

c 

E 
a 



0.51 

8 

1.38 

1.64 

72 

0.17 

1.45 



10- VK 

10- Vk 

W/mK 

GJ/Km3 

GPa 



These values are taken from [9], with the exception of Prj,^ q^ noted above. 



CONCLUSION 



Thermo-optic noise results from thermal fluctuations in the coatings used to make high-reflection mirrors. These 
thermal fluctuations affect the measured position of a mirror through the thermo-elastic and thermo-refractive mech- 
anisms. While both of these mechanisms have been known for some years, they were not treated coherently. The 
coherent treatment presented herein shows that these two mechanisms appear with a relative negative sign. The 
effect is to essentially reduce thermo-optic noise to the point of insignificance for second generation gravitational-wave 
antennae. While it is true that our current knowledge of the properties of coating materials is imprecise, it seems 
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unlikely that better measurements, while desirable, will bring thermo-optic noise back into the realm of relevance. 
This fact should help to guide coating research in the coming years. 

APPENDIX A: COATING AVERAGE PROPERTIES 

Optical coatings are made from alternating layers of materials with different refractive indices. For properties other 
than the refractive index, as long as the length scales involved {rx and rg) are large compared to the layer thickness 
(typically < A/2), we can use suitably averaged material properties to represent the coating. The equations given in 
this section are all taken from [9], and are repeated here only for completeness and clarity. 

The thermal expansion coefficient for a given layer k in the coating is 



1 + o-s 
ttfc = Cfk 

1 - CTfc 



— + (1 - ^cTsJ^r 

1 + o-s Es 



(Al) 



and the volume average coefficient for a coating with N layers each of thickness dk is 



fc=l 



where d is the total coating thickness 

N 



d = ^4. (A3) 



fc=i 

To compute the correction factor in section III the average thermal properties of the coating are needed. The heat 
capacity is a simple volume average, 

N , 

k=l 

while the average thermal conductivity involves the inverse 

APPENDIX B: REFLECTION PHASE OF A MULTI-LAYER COATING 

In this appendix we describe the method we use for computing the reflection phase of a multi-layer coating. 
The initial discussion is somewhat pedantic, but it serves to give us a consistent notation which we develop in the 
subsections detailing the thermo-elastic and thermo-refractive mechanisms. 

We start with the effective reflectivity of the interface between materials with refractive indices ni and 712, passing 
from material 1 to material 2, 

ri,2 = — - — . (Bl) 

ni + 712 

Given two such transitions, from 1 to 2 and from 2 to 3, can equate the reflectivity to that of a two mirror cavity 

-r2,i + r2,3 e-"^^ 
''1.2,3 = -j TT^ B2) 

where (f>2 is the round-trip phase in material 2. 

Note that the reflectivity of the 2 to 1 transition appears in (B2) with the indices in the order seen from inside the 
cavity. In the following text we will use the relation rfe+i.fc = —rk^k+i, to keep the indices in increasing order, and 
then drop the second index, such that r^ = ru^k+i- 
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If we number the interfaces in our coating in the order of increasing depth (i.e., the coating layer in contact with 
the vacuum is 1, and the layer in contact with the substrate is N) we can define a recursion relation using (B2) 

1 + rfcffe+i ^ ' 

where r^ = e~^'^''rk,k+i,...,N is the effective reflectivity of a coating layer, including the round-trip in that layer. The 
base case for this recursion relation is the transition from the A''*'' coating layer to the substrate, 

f^ = e-"f'«rN,s (B4) 

which can be evaluated with (Bl). 

Extending our coating to include the external vacuum as layer provides a natural end to the recursion. The 
reflectivity of the coating is then given by re = fo, and we can use 0o to account for the overall expansion of the 
coating into the vacuum with 0o = Ac/dl, where Ac is the total change in coating thickness. 

To use (B3) to compute changes in reflection phase one must take the derivative with respect to the round-trip 
phase in each layer. Here we give the recursion relation and base case for these derivatives, 

^ n-l-r,.r,_ , .y-! rtih, 1^ ^ J 



d^,-^-'''^ fc=.7 (B5) 



H 

,0 k>i 

From the derivatives of the reflectivity of each layer, the derivative of the reflection phase of the coating as a whole is 

a^ ^ 9arg(fo) ^ -j-f}_^\ Cgg^ 

dcjjk d(l)k \ i~o d(l)k ) ' 

For any quarter or half-wave coating, tq is entirely real and its phase derivatives are entirely imaginary, so much of 
the apparent complexity is not real. 

1. Thermally Induced Changes 

For phase changes induced by a uniform change in temperature we have 

dT ^ d4>k dT' ^ ' 

k=0 ^ 

The phase change due to thermo-elastic and thermo-refractive effects in a coating with layers of thickness dk are 

-^ = ^(Pfc +aknk)dk = -j-Bkdk 

"^ ^ fe=i "'t' 

where, as previously noted, we use 0o to account for the overall expansion of the coating. For any real coating, one 
can evaluate this expression numerically, and thus find d(j)c/dT for that coating. 

2. Relative Sign of TE and TR in 1/4-wave Coatings 

Of particular interest are high-reflection coatings made of 1/4-wave layers of alternating low-n and high-n material. 
For simplicity, we'll assume that the high-n layers have nn > rig and that the low-n layers have n^ ~ ng. Thus, the 
reflectivity from high-n to low-n, is 

TH = ■ . (B9) 

riL + nH 
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As a transition from the vacuum, the first layer is of low-n material and 1/2- wave in optical thickness, such that 



-^ItS^ <«'»' 



Summarizing, these coatings have the following properties 



1 fc <= 1 
-1 fc> 1 



ro k = 
i^k = { rn k even 
-rn k odd 



From the above we can at least determine the signs of the various phase derivatives. We start by noting that 

sign(ffe) = 




and that (B5) inverts the sign of the derivative for each layer with fc > 1. Even numbered layers start with dfk/d4>k 
positive, experience fc — 2 sign inversions, and thus end with a positive sign. Odd numbered layers, on the other hand, 
start with dfk/d(f>k negative, experience an odd number of sign inversions, and thus these also end with a positive 
sign. Since fo is negative, we are ensured that 

sign (^] = -1, for all fc. (BU) 

\0(pkj 

It follows that, for any high-reflection coating of this construction, thermo-elastic and thermo-refractive effects will 
appear with opposite sign in (B7), thanks to the relative minus sign in (B8). 

3. Approximation for High Reflectors 

While equations B7 and B8 are accurate and easy to use in numerical computation, they offer little intuitive 
understanding and fail to provide a concise expression for the thermo-optic mechanisms. To address this, we give an 
approximation which is useful for high-reflection coatings. 

The thermo-elastic mechanism, which arises from motion of the coating's surface, is accounted for by the fc* = 
term in equation B7 (also the second line in equation B8). This term can be expressed in terms of the average coating 
expansion coefficient as 

dzTE c.z'^4'0 - , m1o^ 

„^ = ol—— — acd B12 

The thermo-refractive mechanism is accounted for by the terms with fc > in (B7), which can be thought of as the 
change in reflection phase as measured at a point on the coating's surface. We define an effective TR coefficient /3 
such that [27] 

^ = -/3A (B13) 

For a coating made entirely of 1/4-wave doublets, (3 can be approximated by 



n\l3H + n]jl3L 
i{nl - nl) 



Pq^^ '-^^^Y "27 ^ (B14) 



as given in [6]. A 1/4-wave cap layer is, however, counter-productive and not used in high-reflectors. 
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To find /3 for the common HR coating (made of 1/4- wave doublets with a 1/2- wave cap layer), we modify /3qw by 
approximating f^ ~ sign(ffc). Since the sign of fp is minus in the 1/4-wave case and plus in the 1/2-wave case, each 
term in (B5) with j > is reduced by 



{l + ror _ 1 



n A2 - 2 ■ (B15) 

(1 - roY nj 



Furthermore, we must include the additional 1/4-wave of material in the thicker cap layer 

d(l>c,HW l + ra 1 



so that the additional temperature sensitivity is 



H 



lJc,HW 



1 - ro UL 



1 ^ 1 



(B16) 



Putting these corrections together gives 



^^^^^r'^^,- ^""''^ 



n| 4ni 



which can be rearranged to match equation 5. Equation 7 arises simply from the sum of TE and TR terms 

For alternating layers of Si02 and Ta205, this approximation is within a few percent for coatings with more than ~ 6 
doublets. 

APPENDIX C: RELATIONSHIP TO SUBSTRATE THERMO-ELASTIC NOISE 

The spectrum of thermal fluctuations described by (1), and derived previously in [7] and [8], can be rearranged 
with the help of the thermal diffusion length. If we rewrite (1) as 



S. 






we can see that the first fraction is the spectral density of the thermodynamic fluctuation in a volume defined by the 
diffusion length, while the second is the Gaussian beam average over these volumes. 

The coating thermo-optic coupling is designed such that a similar equation applied to the substrate would result 
in zero. The reason for this is that the loss associated with the coating results from non-adiabatic heat flow due to 
the difference between the coating and substrate. The substrate thermo-elastic noise, on the other hand, results from 
adiabatic heat flow on the scale of the the beam radius re, and is thus smaller by a factor of ~ rr/rc- See, for 
instance, equation 2 in [7] which can be written in our notation as 

= S^l ^/2^!l(a,rT)2 (CI) 

To give an idea of the relative importance of substrate and coating thermo-elastic noise, we divide the coating 
thermo-elastic noise in (26) by (CI) and define the thermo-elastic ratio 

In the case of a gravitational- wave interferometers we have roughly, rg ~ 5 cm and d ~ 5 /xm. For a fused silica 
substrate, rx ^ 40 /im around 100 Hz, such that Rte ~ 10(Aa/Q;s)^ ~ 150, indicating that the substrate contribution 
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is insignificant. For a sapphire substrate [rx ~ 130 /im and as = 5.6 x 10~^), on the other hand, the substrate 
contribution is dominant Rte ^ 0.2{Aa/as)'^ ~ 0.1. 
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